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EFFECT OF A GAS-SOLID SUSPENSION WORKING FLUID ON 
RADIATOR CHARACTERISTICS FOR A 1-MEGAWATT 
BRAYTON SPACE POWER CYCLE 
by Robert Pfeffer,* Salvatore Rossetti, t and Seymour Lieblein 

Lewis Research Center 

SUMMARY 

An analysis has been conducted of the radiator heat-transfer and geometric charac- 
teristics of a high-temperature 1-megawatt single-loop Brayton space power cycle con- 
taining inert-gas graphite -Suspension' working fluids. The analysis considered helium, 
neon, and argon with inclusion of internal radiation heat transfer from the suspension 
particles. Cycle temperature ratios and particle loading ratios for minimum radiator 
planform area were determined for a tube-flow Reynolds number of 50 000. 

For optimum conditions (minimum planform area), -particle radiation had a rela- 
tively small effect on overall heat transfer, and resulted in only a 2 to 3 percent reduc- 
tion in minimum radiator planform area. Maximum indicated reduction in minimum 
planform area compared with the pure gas case ranged from 9 percent for the argon sus- 
pension to 13 percent for the helium suspension for the range of variables and cycle con- 
ditions covered. Considerably larger relative variations were found for tube number, * 
length, and internal diameter and for panel aspect ratio. Tube and panel geometry also 
varied substantially with the gas used. 

The overall reduction in planform area resulting from the use of the gas- solid sus- 
pension varied with the magnitude of the cycle component efficiencies. For identical 
component efficiencies, the planform area reduction calculated for the high-temperature 
1-megawatt cycle was nearly the same as that for a comparable low-temperature 
100-kilowatt cycle analyzed in an earlier study. 


Associate professor, City College of City University of New York. 
^Graduate student, City College of City University of New York. 



INTRODUCTION 


Advanced space power generation systems muSi oe able to continuously supply elec- 
trical power for long periods of time. One such power generating system currently 
under technical evaluation is the indirect-conversion closed- loop heat engine based on the 
Brayton cycle. In this system, heat is generated in a nuclear or solar source and re- 
jected by a radiator, with power being obtained from a turbine located in the working 
fluid cycle. 

Several analyses of potential performance of a Brayton cycle using an inert gas as 
the working fluid have been conducted (refs. 1 to 4). These studies have indicated that, 
due to the relatively low temperature levels in the radiator and the inherently low heat- 
transfer coefficients of gases, the Brayton- cycle waste-heat rejection structure in a 
single- loop system is relatively large and heavy. 

One method that has been proposed to decrease the radiator area and, hence, its 
weight in a single-loop system is the use of a gas-solid suspension rather than a pure 
gas as the cycle working fluid. The effect of using a gas-solid suspension as the working 
fluid on the radiator characteristics of a 100-kilowatt Brayton space power generation 
system has shown (ref. 5) that the addition of graphite particles to an inert gas working 
fluid can result in a reduction of approximately 10 to 20 percent in the required minimum 
radiator planform area. 

As the design power level of the Brayton cycle powerplant is increased, the average 
temperature in the radiator generally increases. Therefore, at high power levels, 
direct radiation from the solid particles to the tube wall (which was neglected for the 
100-kilowatt system studied in ref, 5) might conceivably increase the overall heat trans- 
fer to the wall and thus further lower the required radiator planform area. 

This study determines the overall effect of gas- solid suspensions and particle radi- 
ation on the geometric characteristics of the radiator of a 1-megawatt Brayton space 
power cycle. The analysis is based on a radiator model identical to that used in refer- 
ence 5, but with the addition of particle radiation to the tube wall in the heat-transfer re- 
lations. The effect of particle loading ratio on the convection and radiation heat- transfer 
coefficients and on the radiator planform area and panel tube geometry is also presented. 


ANALYSIS 

Geometrically, the radiator can be identified by its planform area App tube inside 
diameter Dp tube length L, number of tubes N, and panel aspect ratio <p (symbols are 
defined in appendix A). These radiator geometric characteristics are largely determined 
•by the thermodynamics of the power cycle and the heat-transfer and pressure drop capa- 
bilities of the working fluid. The governing equations that are required in calculating the 
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effect of various gas- solid suspension working fluids on these geometric characteristics 
of the radiator are presented herein. The procedure used in developing these relations 
is basically the same as that described in detail in reference 5. 


Power Cycle 

A schematic flow diagram of a single- loop Brayton cycle is shown in figure 1, The 
numbers correspond to the state-point designations used in this analysis. The gas leav- 
ing the heat source at point 1 expands through the turbine to point 2 thereby producing the 
mechanical work necessary to drive the compressor and alternator. From the turbine 
the gas enters the recuperator where it is cooled to point 3, as it transfers heat to the 
gas leaving the compressor. Final cooling to point 4 takes place in the radiator where 
the excess heat is rejected to space. The gas leaving the radiator is then compressed to 
point 5, heated in the recuperator to point 6, and further heated back to point 1 in the heat 
source. 

In general, a power cycle analysis provides equations that enable the determination 
of the working fluid flow rate W, the prime radiating area Ap, and the radiator terminal 
temperatures Tg and T^ from overall cycle parameters and component performance 
• inputs. Prime radiating area Ap is defined as the area of either a tubular radiator 
without fins or a tube-and-fin radiator with 100 percent efficiency. Equations for these 
quantities for the Brayton cycle shown in figure 1 have already been developed in refer- 



Figure 1. - Flow diagram of single-loop Brayton cycle, (Numbers 
are state-point designations.) 
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ences 1 and 5 for the case of no particle radiation heat transfer. However, in this study 
suspension particle radiation to the tube inside walls will be included in addition to heat 
convection. 

When solid particles are present in a flowing gaseous stream which is at a higher 
temperature than the surrounding wall, heat will be transferred to the wall by radiation 
from the particles in addition to heat transfer by convection from the suspension. (It is 
assumed that heat conduction through the tube walls has a negligible effect on the radiator 
external temperatures. ) The radiant heat flux to the wall due to the particle laden flow 
will depend on the gas and wall temperatures and on overall view. factor or radiation effi- 
ciency f , defined as the ratio of the actual radiation emitted by the suspension compared 
with the maximum radiative energy that can be emitted. 

Thus, for an element of tube length, the heat transferred from the gas- solid suspen- 
sion by both convection and radiation to the tube wall must equal the heat radiated to 
outer space, so that 

h i( T s • T wK + - T wK = «( T i - T sbK (1) 

The determination of the value for & is discussed in appendix B. 

The convection heat-transfer coefficient hj^ based on prime radiating area is 
defined as 


h 


R 



( 2 ) 


where for a suspension 


h i 



Similarly, the radiation heat-transfer coefficient h^ based on prime radiating area Ap 
is defined as 


dA. 

h Q = fo — L 
“ “Ap 

Substitution of equations (2) and (3) into equation (1) yields 


( 3 ) 


h R ( T s - T > h q( T s -<) = «(' 


IT’ - T 4 ) 
1 w sp/ 


(4) 


Equation (4) relates the wall temperature T w to the suspension bulk temperature 



T g at any point along the radiator. For a given radiator outlet temperature T^, the 
standard cycle analysis (ref. 1) determines the radiator inlet temperature Tg. The 
radiator temperature drop (Tg - T^) can then be arbitrarily divided into small segments 
and the wall temperature at each of these segments can be determined by solving equa- 
tion (4) numerically. Once the wall temperatures have been determined, the prime radi- 
ating area A can be determined from the following heat balance around an element of 

ir 

radiator tube length: 


so that 



(5) 


( 6 ) 


Equation (6) can be integrated numerically along the points where the wall temperature 
has been computed between the inlet and outlet temperatures of the radiator. 


Gas-Solid Suspension 

The gas- solid suspension characteristics required in this analysis were taken mostly 
from the relations presented in reference 6. Some of the more important of these rela- 
tions are summarized in the following paragraphs. 

Specific heat and its ratios . - The specific heat Cp of a gas is a unique value for 
that gas. However, for a gas-solid suspension, the value of the specific heat will depend 
on the specific heats of the carrier gas and the particle and on the particle loading ratio 
77 defined as the ratio of the mass flow rate of solid particles to that of the pure gas, 

W„ 

V=-^ (7) 

W„ 
g 

The specific heat of the suspension is given by 

( c pH c p) g (nf) (8) 

where 6 is defined as the ratio of specific heat of the particles to that of the pure gas 


5 



(9) 


6 = 



The proper value to be used for the isentropic heat ratio of a suspension y is not as 

s 

clearly defined as for the specific heat. If, however, equilibrium between gas and parti- 
cles can be assumed (T = T , V = V ), the isentropic specific heat ratio of a suspen- 

& A & 

sion y g in terms of the readily available properties of the system can be estimated by 
(ref. 6) 


/ 1 + 5n 

y g = >v — 

s V + yJv / 


do) 


For a suspension, the specific heat ratio will be less than that of the carrier gas. For 
large loading ratios 77 > 10, y c — 1.0, and the suspension behavior approaches that 
of an incompressible fluid. 

Friction factor. - The friction factor ratio fl/f„ is assumed to be adequately repre- 

S g 

sented by the relation given in reference 6 


f’ 

- = 1 + 4. 0(Re) 

f ° 

g 


-0.32 


( 11 ) 


For fully developed turbulent flow in smooth pipes 


f = 0.046 (Re) 

o 


- 0.2 


( 12 ) 


Convective heat transfer. - The convective heat- transfer coefficient ratio h c /h„ 

to g 


based on tube inside area is given by (ref. 6) 

h 
h 


— = 1 + 4. 0(Re)" 0, 32 677 

O 


(13) 


g 


For turbulent flow 


h =Q.023 -£ (Re)°* 8 (Pr)°* 3 

& D. & & 


or 
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h = 0. 023 -£ (Re)?’ 8 


D i 


( c p). 


M c 


£. 


0.3 


(14) 


Suspension radiative heat transfer . - The suspension radiation efficiency in equa- 
tion (1) is determined from the equation (appendix B) 


i = J_ T + _L-I + JL 

^ 16 e w 2 8 t 


(15) 


where the optical thickness r is approximated by the relation 


3 Pg D i 

T = - 7] — 

2 p p d p 


(16) 


In order to compare the effect of suspension radiation heat transfer with that of the 
convection heat transfer, a pseudo-radiation heat-transfer coefficient hQ is defined as 


h Q 



(17) 


Since T g and T w vary along the radiator tube, an average value of hQ was deter- 
mined by applying equation (17) to a number of equally spaced points along the tube and 
averaging the values obtained. The relative magnitude of the effect of radiation from the 
particles to the tube wall compared with the total heat-transfer rate to the tubes is indi- 
cated by the quantity x}/ defined as 


* = h Q- (18) 

h R + h Q 

Radiator Geometry 

The radiator configuration considered in the analysis is identical to that presented in 
reference 5 and consists of a central fin-tube panel geometry of axial length L and N 
number of tubes as indicated in figure 2. The tubes are of inside diameter and out- 
side diameter D Q , and the rectangular fins are of thickness tj. and width 2 Heat re- 
jection occurs from the outside surface of the tubes and from both surfaces of the fins by 
radiation to space. For simplicity, the headers are neglected, and only the fin-tube panel 
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Tube 



(b) Cross sectional view. 

Figure 2. - Central fin -tube radiator configuration. 


is considered. In addition, representative constant values of l^/ D Q , t^, and D q /D^ are 
chosen for all cycle conditions as was done in reference 5. Based on this simplified 
model geometry and these assumptions, equations determining the radiator pressure re- 
lations and radiator panel geometry were developed in reference 5. The principal equa- 
tions will be repeated in the following paragraphs. 

Heat-transfer coefficient. - The cycle calculation requires a value of convective 
heat-transfer coefficient h R’ based on prime radiating area A p , given by equation (2). 

In terms of the radiator model of figure 2, the prime and inside areas and A i of 
.equation (2) are, respectively, 


A 


P 



1 + 2rjj 



LN 


(19) 


and 


A i = ttDjLN 


( 20 ) 
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Thus, the convective heat-transfer coefficient h^ is (from eqs„ (2), (19), and (20)) 




( 21 ) 


For a suspension h. = h_, where h_ is obtained from equations (13) and (14). 

IS s 

Tube geometry . - The tube inside diameter is obtained from the relation for 
friction pressure drop (ref. 5) based on equations (11) and (12), and the relation for gas- 
flow Reynolds number given by 



( 22 ) 


so that 




■ [l + 4. 0(Re)” 0 ' 32 
( AP >rad L g 


V 


-|/LFjU 2 R'T 3 




c 3 


where 



and 


P 3 




( AP >rad 


(23) 


(24) 


(25) 


In equation (25), K 2 is defined as the ratio of the pressure drop in the cold side of the 
recuperator to the total pressure drop in the cold side of the cycle (radiator and recuper- 
ator), or 

( AP ) rec , „ 

K„ = (26) 

2 ( AP ) C 
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Furthermore, the value of P2 is obtained from the turbine pressure ratio given by 



( 27 ) 


For equations (23) and (25), the pressure drop across the radiator 
tained in terms of the cycle loss pressure ratio from 


<aP >rad 18 ob - 


< iP W = 



(28) 


where is defined as the ratio of the pressure drop in the radiator to the total system 
pressure drop, or 


K l s 


(AP) 


rad 


(AP) 


tot 


(29) 


The tube outer diameter D Q is then 


D = D. 
o 1 



and the tube length L is obtained from equations (19) and (22) as 

A-CReL/oKl + r?) 


L = 


'g' g 




(30) 


(31) 


Panel geometry . - The number of tubes in the panel N is obtained from equa- 
tion (22) as 


N = 


4W s 

^ g (Re) g D.(l + 17) 


(32) 


and the fin length is 
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(33) 


l 


f 



The planform area of the radiator A pl is defined as the projected area of the radiator 
panel 



(34) 


The panel aspect ratio cp is defined as the ratio of overall width to length of the radi- 
ator, or 


(35) 

N(D 0 + 2* f ) A pl 

CALCULATIONS 
Method of Solution 

Once the cycle parameters 77 ^, 77 ^, E, and r^/r^, were specified, and the inde- 
pendent variables Tg/T-p T^/T^, P^, and P g were chosen, the cycle was fully de- 
scribed for a given pure gas or a specified gas-solid suspension. From the cycle anal- 
ysis, W and all gas temperatures were determined according to the corresponding equa- 
tions given in reference 5. Suspension radiation efficiency P was obtained using equa- 
tion (15). The wall temperatures along the radiator tubes were computed at 21 equally 
spaced points by solving equation (4) using the Newton- Raphs on technique. Prime radi- 
ating areas A were obtained by integrating equation (6) numerically using Simpson’s 

ir 

method. 

By fixing the radiator model constants Dq/Dj, Z f /D Q , 77 ^, Kj, and Kg and by pre- 
scribing the gas Reynolds number (Re)^ and loading ratio 77 in the tubes of the radiator, 
equations (23), (31), (32), (34), and (35) in conjunction with the cycle outputs were used 
to determine D^, L, N, A^, and cp, respectively. An iterative procedure was required 
based on an initial assumption of D i because h R (which is necessary for the cycle anal- 
ysis) depends on 

The minimum value of the planform area A^ rather than the prime radiating area 
Ap was selected as the criterion for determining both the optimum cycle temperature 
ratios and the optimum suspension loading ratio tj for each gas. The planform area was 
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considered to be more representative of radiator size and weight than prime area. Opti- 
mum temperature ratios were determined from a scanning of the computer outputs for 
each value of 77. 


Inputs 


The example chosen for the parametric analysis was a 1-megawatt shaft-power out- 
put cycle with a suspension of graphite in either pure helium, neon, or argon as the 
working fluid. From turbomachinery considerations, compatible values of = 3000° R 
(1667 K) and P ^ = 5. 8752x10^ pounds per square foot (2. 813 MN/m^) were chosen as the 
turbine entering temperature and pressure, respectively. Values of the compressor, 
turbine, and recuperator efficiencies were taken as 0. 83, 0. 89, and 0. 80, respectively, 
for both the pure gases and the gas- solid suspensions. The independent cycle parameters 
are Tg/Tj, T^/T^, and the loss pressure ratio r T /r^. In this analysis, T 2 /T^ was 
varied between 0. 70 and 0. 95, and T^/T-^ was varied between 0. 25 and 0. 50, both in in- 
crements of 0.025. All the calculations were conducted at r^/r^ = 0.90, Kj = 0.5, and 
K 2 = 0. 25. 


The properties of the inert gases (helium, neon, and argon) required as inputs to the 
analysis were the thermal conductivity k , the viscosity u , the specific heat (C n ) , 

o o Jr g 


the isentropic specific heat ratio y , and the gas constant R\ Values of these proper- 
ties for the temperature range of interest are presented in appendix C. The gas-solid 
suspensions studied were suspensions of 30-micrometer-diameter graphite particles in 
these same inert gases. The pertinent physical properties of graphite are also presented 
in appendix C. Loading ratios 77 were varied from 0 to 7 in increments of 1 for helium- 
graphite, from 0 to 3 in increments of 0. 5 for neon- graphite, and from 0 to 1. 5 in incre- 
ments of 0.25 for argon- graphite. 

Based on previous analyses of gas cycle radiators, representative geometric param- 
eter values were selected as D q /D^ = 1. 33, £j/D 0 = 3.5, and 77^ = 0. 67. For simplicity, 
a single gas Reynolds number of 50 000 was selected for the calculations. The equivalent 

sink temperature of the radiator T_._ was chosen as 400° R (222 K), and a value of 0. 86 

°P 

was chosen for the emissivity e of both the inner and outer surfaces of the tubes and of 
the fins. 


RESULTS AND DISCUSSION 
Optimum Temperature Ratios 

Table I summarizes the principal calculated results of the analysis for the three sus- 
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TABLE I. - RADIATOR CHARACTERISTICS FOR 1-MEGAWATT BRAYTON CYCLE WITH PARTICLE-TO-WALL RADIATION 


[Optimum values at gas Reynolds number of 50 000.] 


Loading 

ratio, 

V 

Optimum temperature 
ratios 

Radiator 
planform 
area, A pl 

Convective heat- 
transfer coef- 
ficient, h R 

Pseudoradiation 
heat-transfer 
coefficient, hq 

Tube inside 
diameter, 

Tube length, 
L 

Number 

of 

tubes, 

N 



Panel 

aspect 

ratio, 

<P 

Turbine 
inlet to 
outlet, 
V t i 

L 

Compressor 
inlet to 
turbine 
outlet, 

V T i 

ft 

m 

ft 

m 

ft 2 

rvr 2 

m 

Btu 

W 

Btu 

W 

(ft 2 )(°R) 

<m 2 )(K) 

(ft 2 )(°R) 

(m 2 )(K) 

Helium-graphite 

0 

0.775 

0.350 

2626 

244 

58.2 

330.3 

0 

0 

0.071 

0.022 

54.4 

16.6 

63.9 

1.13 

1 

.800 

.375 

2427 

225 

53.3 

302.4 

1.40 

7.94 . 

.084 

.026 

58.5 

17.8 

46.4 

1.41 

2 

.800 

.375 

2337 

217 

42.9 

243.4 

2.08 

11.80 

.107 

.033 

74.8 

22.8 

27.4 

2.39 

3 

.825 

.400 

2294 

213 

42.3 

240.0 

2.63 

14.92 

.116 

.035 

73.7 

22.5 

25.3 

2.37 

4 

.825 

.400 

2281 

212 

35.7 

202.6 

2.68 

15.21 

.142 

.043 

88.5 

27.0 

17.1 

3.43 

5 

.825 

.400 

2301 

214 

30.4 

172.5 

2.61 

14.81 

.173 

.053 

104.7 

31.9 

12.0 

4.75 

6 

.850 

.425 

2313 

215 

33.0 

187.3 

2.80 

15.89 

.168 

.051 

95.8 

29.2 

13.5 

3.97 

7 

.850 

.425 

2342 

218 

28.9 

164.0 

2.65 

15.04 

.199 

.061 

110.0 

33.5 

10.1 

5.17 

Neon-graphite 

0 

0.775 

0.350 

2951 

274 

32.4 

183.9 

0 

0 

0,039 

0.012 

18.9 

5.8 

373.5 

0.121 

.5 

.825 

.375 

2653 

246 

26.1 

148.1 

1.58 

8.97 

.054 

.016 

27.4 

8.4 

166.9 

.284 

.75 

.825 

.400 

2626 

244 

22.5 

127.7 

2.13 

12.09 

.068 

.021 

30.2 

9.2 

119.9 

.347 

1.0 

.850 

.400 

2632 

245 

22.6 

128.2 

2.31 

13.11 

.070 

.021 

33.7 

10.3 

104.0 

.432 

1.5 

.875 

.425 

2710 

252 

21.8 

123.7 

2.69 

15.26 

.081 

.025 

36.9 

11.2 

85.7 

.503 

2.0 

.875 

.425 

2799 

260 

17.3 

98.2 

2.58 

14.64 

.110 

.034 

48.9 

14.9 

49.0 

.855 

2.5 

.900 

.450 

2938 

273 

19.7 

111.8 

2.72 

15.43 

.106 

.032 

46.6 

14.2 

56.1 

.738 

3.0 

.900 

.450 

3048 

283 

16.6 

94.2 

J 

2.52 

14.30 

.135 

.041 

57.6 

17.6 



37.0 



1.09 

Argon-graphite 

0 

0.775 

0.350 

3133 

291 

25.7 

145.8 

0 

0 

0.019 

0.006 

7.55 

2.3 

2036.0 

0.018 

.250 

.825 

.375 

2855 

265 

20.8 

118.0 

.92 

5.22 

.026 

.008 

11.1 

3.4 

911.8 

.043 

.375 

.850 

.400 

2850 

265 

21.0 

119.2 

1.38 

7.83 

.028 

.009 

11.3 

3.4 

846.5 

.045 

.500 

.850 

.400 

2853 

265 

18.0 

102.1 

1.63 

9.25 

.034 

.010 

13.8 

4.2 

569.2 

.067 

.750 

.875 

.425 

2937 

273 

17.5 

99.3 

2.21 

12.54 

.039 

.012 

15.2 

4.6 

470.3 

.079 

1.00 

.875 

.425 

3033 

282 

14.0 

79.4 

2.37 

13.45 

.052 

.016 

20.1 

6.1 

271.1 

.133 

1.25 

.900 

.450 

3154 

293 

16.0 

90.8 

2.78 

15.78 

.050 

.015 

19.0 

5.8 

311.8 

.114 

1.50 

.900 

.450 

3259 

303 

13.6 

77.2 

2.77 

15.72 

.063 

.019 

23.4 

7.1 

207.2 

.168 


pensions considered. The temperature ratios across the turbine Tg/Tj and the ratios 
of compressor inlet to turbine inlet temperature T^/T^ listed in the table are the tem- 
perature ratios at each value of loading ratio 77 which result in a minimum value of the 
planform area A^. As indicated previously, these optimum temperature ratios were 
determined from a scanning of the computer outputs for each value of 77. 

As can be seen from table I, both optimum temperature ratios increase with increas- 
ing particle loading ratio r\. The rate of increase is smallest with helium and most pro- 
nounced with argon. The increase in the optimum temperature ratios with increasing 
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loading ratio r\ is due to the fact that the isentropic specific heat ratio y 0 decreases 

s 

with loading ratio so that the temperature drops across both the turbine and compressor 
de crease o The variation of the optimum temperature ratios with loading ratio was the 
same as that for the 100-kilowatt cycle studied in reference 5. 

Radiator outlet temperature at the loading ratio 77 for the smallest planform 
area A pl was 1200° R (667 K) for all gases. Radiator inlet temperature Tg for this 
condition was 1770° R (983 K) for helium and neon, and 1735° R (964 K) for argon. 


Heat-Transfer Characteristics 

Convective heat-transfer coefficient . - The effect of the gas-solid suspension on the 
convective heat-transfer coefficient h R based on prime radiating area A p is illustrated 
in figure 3. The values of h R for the suspension decrease with increasing loading ratio 



Figure 3. - Variation of convective heat -transfer coefficient with weighted 
specific heat ratio. 


rj for all three gases. This apparent anomaly, as discussed in reference 5, is due to the 
relative variations of hj for the pure gas and for the suspension. Although the ratio 
h s /hg increases for a given tube diameter D i as loading ratio tj is increased (as indi- 
cated by eq. (13)), the tube diameter for the suspension Dj also increases compared 
with the pure gas case (table I). There is a net decrease in the value of h g compared 
with hg. And with Aj/A p a constant in all cases (this ratio depends only on geometry), 
there is a corresponding net reduction in h R for the suspension as compared with the 
pure gas case. Figure 3 also shows that the convective heat-transfer coefficient h R is 
highest for helium- graphite suspensions and lowest for argon-graphite suspensions. 

Radiative heat-transfer coefficient . - The heat transfer by direct radiation between 
the particles and the tube wall is dependent on the radiation efficiency 3F and is also a 
function of the difference between the suspension bulk temperature and the radiator wall 
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With radiation 

Without radiation 



(a) Helium-graphite, particle loading ratio, 7.0. 



2 .3 .4 .5 .6 .7 .8 


Fraction of tube length from entrance 
(b) Argon -graphite; particle loading ratio, 1.5. 
Figure 5. - Variation of temperature along tube length. 




temperature (eq, (1)). Since the emissivity of the wall e is assumed constant at 0. 86 , 
& depends only on the optical thickness r, which will change with different cycle condi- 
tions and particle loadings. Figure 4 shows the variation of radiation efficiency with 
particle loading ratio 77 for the helium, neon, and argon suspensions using 30- micro meter- 
diameter graphite particles at optimum cycle conditions (optimum Tg/T^ and T^/T^). 
For all cases considered, the radiation efficiency & rises rapidly as particle loading 
ratio 77 is increased, reaches a maximum value, and slowly decreases. For the wall 
emissivity used in this analysis, equation (16) predicts a maximum view factor 2F of 
0. 634. Figure 4 shows that maximum & is approached at loading ratios 77 of approxi- 
mately 4.0, 1.5, and 1.5, respectively, for the helium-, neon-, and argon-graphite sus- 
pensions. 

Figure 4 also shows the variation of the radiation heat-transfer fraction i// with 
loading ratio 77. Although & tends to reach a maximum with loading ratio 77 , the quan- 
tity 1 }/ continues to increase, but at a reduced rate, for the range of loading ratios 77 
covered for each gas. At the maximum values of F the contribution of the radiation heat 
transfer to the total heat transfer 1 )/ was around 7. 5 percent for helium, 11.5 percent 
for neon, and 17 percent for argon. 

Figure 5 illustrates the effect of particle to wall radiation on the radiator tube wall 
temperature at the largest values of 77 considered for the helium- graphite and argon- 
graphite suspensions. In both cases the difference between the wall temperature, when 
taking into account particle radiation or when neglecting it, is relatively small, varying 
from 20° to 50° R (11 to 28 K) at the radiator inlet and from 3° to 15° R (2 to 8 K) at the 
radiator exit. This small variation in temperature, when compared with the difference in 
temperature between the wall and the bulk suspension (as also shown in fig. 5), is a re- 
flection of the moderate radiation contribution to the total heat-transfer rate, as already 
revealed in figure 4. 


Radiator Geometry 

Tube geometry . - The effect of loading ratio 77 on such panel geometry character- 
istics as the length, number, diameter of radiator tubes and radiator panel aspect ratio 
can be obtained from the values listed in table I. The table shows that, although the plan- 
form area A^ decreases to a minimum and then increases as the loading ratio 77 is in- 
creased, the diameter and length L of the radiator tubes continuously increase, 
whereas the number of radiator tubes N continually decreases in the loading ratio range 
of interest. The variations in D^, L, N, and ip were quite pronounced over the range 
of particle loading ratios covered, with maximum changes in N and cp of around an 
order of magnitude for the neon and argon suspensions. At minimum planform area, the 
tube length, inside diameter, and number varied from 50 to 100 percent of the corre- 
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sponding values for the pure gas, while the panel aspect ratio cp increased by approxi- 
mately a factor of 3. The use of a solid suspension therefore provides a very effective 
means for altering the tube and panel geometry at a given gas-flow Reynolds number. 

The table also shows that, for a fixed value of gas-flow Reynolds number, the use of 
helium results in a panel aspect ratio cp approximately 10 times that of neon and 100 
times that of argon at minimum A^. These same trends were observed for the 100- 
kilowatt cycle studied in reference 5. 

Planform area. - The variation of panel planform area V with particle suspension 
loading ratio r\ including the effect of particle radiation to the wall is shown by the solid 
curves in figure 6 for the three gases considered. Values of the planform area for 
77 = 0 represent the pure gas case. The figure shows that a minimum planform area Ap^ 
is obtained at a specific loading ratio 77 for each gas. The optimum loading ratios r\ 
are approximately 4.0, 0. 85, and 0.375 for the helium-, neon-, and argon- graphite sus- 
pensions, respectively. Thus, only for helium do the radiator planform area Ap^ and 
radiation efficiency & optimize at the same value of loading ratio. 

The dashed curves plotted in figure 6 represent the planform area Apj as a function 
of loading ratio neglecting particle -to- wall radiation. Comparing the dashed and solid 
curves shows that the optimum loading ratio t] is slightly higher when the particle radi- 
ation is considered but results in only 2 to 3 percent decrease in planform area Ap^ at 
the optimum loading ratios rj. The decrease in planform area resulting from particle 
radiation tended to increase for values of rj greater than optimum, especially for the 
neon and argon suspensions. For example, at 77 for maximum the percent reduction 
in planform area due to particle radiation was 4. 5 and 11. 0, respectively, for neon and 
argon. 

A minimum of the planform area Ap^ with loading ratio 77 is basically due to the 
opposing effects of the isentropic specific heat ratio y a and the convective heat-transfer 
coefficient h^. For a given gas, an increase in the solids loading ratio r\ decreases 
y g , which increases the cycle efficiency and thus tends to decrease the planform area 
Apj. At the same time, an increase in 77 decreases h^ (see fig. 3) and, therefore, 
tends to increase the planform area Ap^„ The average radiation heat-transfer coeffi- 
cient hQ, which is directly related to xp , increases with loading ratio 17 for the condi- 
tions considered in the analysis (see fig. 4). However, the contribution of hQ to the 
tube heat-transfer to the wall is only moderate so that the net effect of hQ on the opti- 
mum loading ratio 77 is a slight shift in the direction of higher loading ratios as shown 
in figure 6 . 

Since y g and h-^ are both functions of the weighted specific heat ratio 677 this 
quantity will define the minimum planform area V This is illustrated by figure 7, 
which shows that for all three graphite-gas suspensions considered the minimum plan- 
form area Apj occurs at approximately the same value of 677 equal to 1. 3 at the gas 
Reynolds number of 50 000. Based on the results of figure 7, the magnitude of the de- 
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Weighted specific heat ratio, 6rj 

Figure 7. - Variation of radiator planform area with weighted specific 
heat ratio. 

crease in panel planform area obtainable by the use of a graphite-gas suspension in the 
1-megawatt Brayton cycle is approximately 13 percent for helium, 11 percent for neon, 
and 9 percent for argon; 

In order to compare the reduction in planform area obtained in the 1-megawatt 
cycle with the 100-kilowatt cycle studied in reference 5, the ratio of planform area re- 
quired with a helium -graphite suspension to the planform area required with pure helium 
was plotted against particle loading ratio rj (fig. 8). The upper three curves represent 
the planform area for the 1- megawatt cycle case with and without consideration of parti- 
cle radiation in the tube internal heat transfer. The lowest curve is the variation ob- 
tained for the 100-kilowatt cycle suspension as reported in reference 5. 

It is seen that the percentage reduction in radiator planform area for the 100- 
kilowatt cycle without consideration of particle radiation was actually greater than the 
percentage reduction for the 1-megawatt cycle with radiation included. These results 
were typical of the three suspensions studied. Since the percentage reduction in the con- 
vective heat-transfer coefficient h R was practically the same for both cycles, these 
results can be explained on the basis of the different values of compressor, turbine, and 
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recuperator efficiencies used in the two studies. The differences in component effi- 
ciencies influenced the effect of the isentropic specific heat ratio y e on the overall 
cycle efficiency and the working fluid mass flow rate W , which, in turn, affected the 
prime area of the radiator Ap to produce a smaller reduction in planform area Ap^. 

The variation of planform area ratio for the 1-megawatt cycle with particle radiation 
for the same component efficiency values as for the 100-kilowatt cycle is given by the 
long-dashed curve. A sizable difference is observed compared with the original results. 
The relative effect of the use of a suspension working fluid will therefore depend on the 
magnitude of the component efficiencies selected for the cycle study. 

On the basis of the same component efficiencies, figure 8 indicates planform area 
results that are nearly the same for both cycles considered. Although the calculations of 
reference 5 did not include the effects of suspension particle radiation, the results ob- 
tained can be considered to be very nearly representative of the values if particle radia- 
tion had been included. This judgment is based on the observation of the relatively small 



Figure 8. - Effect of cycle condition on ratio of radiator planform area with sus- 
pension to planform area with pure gas. Helium-graphite system. 


21 



effect of particle radiation on minimum planform area in the high-temperature 
1-megawatt case (2 to 3 percent) and the substantially lower cycle temperatures con- 
sidered for the 100-kilowatt case. Thus, it can be concluded that the planform area re- 
duction results of reference 5 are also representative of the higher power cycle. Maxi- 
mum reduction in minimum planform area found in the study of reference 5 was around 
16 percent for helium. 


SUMMARY OF RESULTS 

The analysis of the application of an inert gas-graphite suspension with particle radi- 
ation as the working fluid in a high-temperature 1-megawatt single- loop Brayton cycle 
has produced the following major results for optimum cycle temperature ratios and a 
tube-flow Reynolds number of 50 000: 

1. Maximum values of particle-to-wall radiation efficiency 3F were attained at par- 
ticle loading ratios t] of approximately 4.0, 1.5, and 1.5, respectively, for the helium, 
neon, and argon suspensions. At these values of radiation efficiency, particle radiation 
contributed approximately 7. 5, 11.5, and 17 percent, respectively, to the total 
suspension-to-tube heat transfer. However, the net effect of particle radiation on mini- 
mum radiator planform area Ap^ was relatively small, resulting in only an approxi- 
mately 2 to 3 percent decrease in radiator planform area compared with the value ob- 
tained when particle radiation was neglected. 

2. The particle loading ratios 17 for minimum radiator planform area Ap^ were 
4.0, 0.75, and 0.375, respectively, for the helium, neon, and argon suspensions. The 
corresponding values of reduction in planform area compared with the pure gas case were 
about 13 percent for helium, 11 percent for neon, and 9 percent for argon. At minimum 
planform area Ap^, the weighted specific heat ratio 677 was the same at 1 . 3 for all three 
gas suspensions. 

3. The use of a gas- graphite suspension had a pronounced affect on the number, 
length, and inside diameter of the tubes and on the panel aspect ratio. Variations of be- 
tween 50 and 100 percent compared with the pure gas cases were indicated for tube length, 
diameter, and number at minimum planform area. Panel aspect ratio was increased by a 
factor of around 3. 

4. The magnitude of the reduction in radiator planform area A p ^ resulting from the 
use of a gas-solid suspension working fluid varied significantly with the magnitude of the 
component efficiencies used in the cycle. 

5. The calculated radiator characteristics and optimum temperature ratio values 
were generally similar to those determined for a comparable low-temperature 100 - 
kilowatt cycle analyzed in a previous study. Furthermore, when identical cycle compo- 
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nent efficiencies were considered, the reduction in radiator planform area A j for the 
high-temperature 1-megawatt cycle was nearly the same as for the low-temperature 100- 
kilowatt cycle. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 17, 1969, 

120-27. 
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APPENDIX A 


SYMBOLS 


A 

o 

particle cross sectional area, ft ; 

2 

m 

h q 

radiation heat-transfer coefficient 
defined by eq. (3), Btu/ 

A i 

radiator tube inside surface area, 


(hr)(ft 2 )(°R) 4 ; W/(m 2 )(K) 4 


ft 2 ; m 2 

h R 

convective heat-transfer coeffi- 

A P 

2 2 

prime radiating area, ft ; m 

2 2 

radiator planform area, ft ; m 

cient based on radiator prime 
area, Btu/(ft 2 )(°R)(hr); 

V 


W/(m 2 )(K) 

C P 

fluid heat capacity, Btu/(lbm)(°R); 
J/(kg)(K) 

K 1 

radiator pressure-drop fraction: 
ratio of pressure drop through 

D i 

radiator tube inside diameter, ft; 
m 


the radiator to total pressure 
drop in the cycle 

D o 

radiator tube outside diameter, 
ft; m 

K 2 

recuperator pressure-drop frac- 
tion: ratio of pressure drop 

D P 

E 

F 

particle diameter, ft; m 
recuperator effectiveness 
correction factor for nonisother- 

k 

through the recuperator to cycle 
cold side pressure drop 

absorption coefficient, 1/ft; 1/m 

9 

h 

mal flow in the radiator 
radiation efficiency or view factor 
Fanning friction factor for a pure 

k g 

L 

thermal conductivity of gas, Btu/ 
(sec)(°R)(ft); J/(sec)(m)(K) 

length of radiator tubes, ft, m 

gas 

*f 

half width of fin, ft; m 


friction factor for gas- solid sus- 

N 

number of radiator tubes 

pension based on gas density, 
defined by eq. (11) 

n 

number of particles per unit vol- 

o q 

ume of suspension, 1/ft ; 1/m 


gravitational constant, (lbm/lbf) 

P 

9 o 

gas pressure, lbf/ft ; N/m 


(ft/sec 2 ); (kg)(m)/(N)(sec 2 ) 

Pr 

Prandtl number 

h i 

convective heat-transfer coeffi- 
cient based on tube inside area, 

P s 

output shaft power, kW 


Btu/ (ft 2 )(°R) (hr); W/(m 2 )(K) 

(AP) C 

sum of pressure drops across ra- 
diator and recuperator (cold 

h Q 

pseudoradiation heat-transfer 
coefficient defined by eq. (17), 
Btu/ (hr ) (ft 2 ) (°R); W/(m 2 )(K) 


pressure drop), lbf/ft 2 ; N/m 2 
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< iP >rad 

pressure-drop across the radi- 

^C 

compressor efficiency 

ator, lbf/ft 2 ; N/m 2 

V f 

fin efficiency 

< aP Vee 

pressure-drop across cold side 
of the recuperator, lbf/ft ; 

7] T 

turbine efficiency 


N/m 2 

P 

gas viscosity, lbm/(ft)(sec); 

< iP >tot 


(N) (sec)/ (m 2 ) 

total pressure- drop through 


density, lbm/ft 2 ; kg/m 2 

heat-transfer components in 

P 


the cycle, lbf/ft 2 ; N/m 2 

p p 

bulk density of particles in sus- 

R' 

gas constant, (lbf) (ft)/ 

pension, lbm/ft 2 ; kg/m 2 


(lbm)(°R); J/(kg)(K) 

CT 

Stefan- Boltzman constant, Btu/ 

(Re) g 

gas Reynolds number 


(ft 2 ) (°R 4 ); J/(m 2 )(K 4 ) 

r T / r c 

cycle loss pressure ratio, 

7 

optical thickness 



<P 

radiator panel aspect ratio 

T 

absolute temperature, °R; K 


particle-radiation heat-transfer 

T sp 

equivalent sink temperature, 


fraction (eq. (18)) 


°R; K 

Subscripts: 

t f 

thickness of fin, ft; m 

g 

gas 

*t 

thickness of radiator tube, ft; 



m 

P 

particle 

V 

velocity through radiator tubes, 

s 

suspension 


ft/sec; m/sec 

w 

conditions at tube wall 

w 

mass flow rate, lbm/sec; 

1 

turbine inlet 


kg/sec 

2 

turbine outlet 

y 

isentropic specific heat ratio 

3 

radiator inlet 

5 

ratio of particle to gas specific 

4 

radiator outlet 


heats 

5 

compressor outlet 

e 

emissivity 

6 

heat- source inlet 

V 

loading ratio of suspension, 




lbm soUd/ lb "W ^soUd/ 
b %as 
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APPENDIX B 


GAS RADIATION EFFICIENCY 

The gas radiation, efficiency or view factor 2F for use in equation (1) was deter- 
mined from Deissler’s diffusion approximation for thermal radiation in gases with jump 
boundary conditions (ref. 7). For a gray gas flowing through a circular tube of constant 
diameter, Deissler’s results may be approximated by the equation (written in the nota- 
tion of this report) 


A. = _L T + _L.I + jl 

P 16 € w 2 8t 


(Bl) 


where r is the optical thickness of the gas and e w is the emissivity of the inside wall 
of the tube. 

In general, for an absorbing, nonscattering, gray gas, 

r = kD. (B2) 


where k is absorption coefficient and D i is tube inside diameter. For particle suspen- 
sions, k is defined as 

k = nA (B3) 

where n is the number of particles in a unit volume of a suspension and A is the aver- 
age radiation absorption cross section of a single particle in the suspension. For par- 
ticle sizes greater than the wavelength of the radiation, it can usually be assumed that 
the absorption cross section of the particle A is equal to its geometric cross section. 
The particle size considered in this analysis (Dp = 30 pm) is sufficiently large to be 
treated in this manner. 

For spherical particles, the number of particles in a unit volume of suspension n 
is given by 


n = 



p ttD 3 

'p p 


(B4) 


where p' p is the bulk density of the particles (lbm solid /ft 3 uspension or kg golid / 

^suspension) and p p is the actual densit y of the P article (lbm solid /ft solid or **aoUd/ 
m solid)‘ However, according to reference 6, the bulk density of the particles Pp may 
be expressed in terms of the loading ratio r\ and the density of the carrier gas Pg by 
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(B5) 


Finally, substitution of equations (B3) to (B5) into equation (B2) results in the relation of 
the optical thickness of a suspension written as 



P P °P 


(B6) 


Equation (B6) together with equation (Bl) was used to determine the gas radiation effi- 
ciency 2F in terms of the inputs considered in this study. 
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APPENDIX C 


PROPERTIES 

The values of the gas properties utilized as inputs to the computer program are pre- 
sented herein. Since both viscosity and thermal conductivity of the gases are dependent 
on temperature, equations to determine these properties as function of the arithmetic 
mean temperature in the radiator are given. These equations were obtained from faired 
curves based on the data of reference 8 in the temperature range from 760° to 1960° R 
(422 to 1090 K). 


Property 

U.S. Customary Units 

SI Units 

Helium 

S 

1.25 Btu/(lbm)(°R) 

5230 J/(kg)(K) 

y 

1.67 

1.67 

r t 

386.25 (lbf/lbm)(ft/°R) 

2080 J/(kg)(K) 

k g 

0.0433 T°- 650 x 10" 5 Btu/(sec)(ft)(°R) 

0. 3950 T 0, 650 xl0 -2 J/(sec)(m)(K) 

M 

0.0238 T 0,650 x 10' 5 lbm/(ft)(sec) 

0.0519 T°- 650 xl0- 5 (N)(sec)/m 2 

Neon 

c p 

0.25 Btu/(lbm)(°R) 

1046 J/(kg)(K) 

r 

1.67 

1.67 | 

R’ 

77.25 (lbf/lbm)(ft/°R) 

416 J/(kg)(K) 

k g 

0.0128 T°- 662 xl0' 5 Btu/(sec)(ft)(°R) 

0. 1176 T°- 662 xl0” 2 J/(sec)(m)(K) 


0.0232 t 0, 665 xl0" 5 lbm/(ft)(sec) 

0.0509 T 0,665 x 10* 5 (N)(sec)/m 2 

Argon 

c p 

0. 125 Btu/(lbm)(°R) 

523 J/(kg)(K) 

y 

1.67 

1.67 

R r 

38.63 (lbf/lbm)(ft/°R) 

208 J/(kg)(K) 

k g 

0.00298 T 0 ' 725 x 10' 5 Btu/(sec)(ft)(°R) 

0.0284 T°- 725 xl0- 2 J/(sec)(m)(K) 


0.0156 T°- 725 x 10‘ 5 lbm/(ft)(sec) 

0.0355 T°' 725 x 10‘ 5 (N)(sec)/m 2 
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The specific heat of graphite over the temperature range 500° to 2000° R (278 to 
1112 K) was taken as 0.400 Btu per pound mass per °R (1674 J/(kg)(K)). Therefore, the 
values of 6 which were used for the suspensions are 


Suspension 

Particle to gas 
specific heat ratio, 


a 

Helium- graphite 

0.312 

Neon-graphite 

1.56 

Argon- graphite 

3.12 


It was also assumed that the average diameter of the graphite particles is 30 micro- 

Q 

meters, the density of the graphite particles is 96.5 pounds per cubic foot (1546 kg/m ), 
and the absorptivity is 1. An average particle diameter of 30 micrometers was selected 
to provide for a stable suspension for the range of optimum loading ratios covered. This 
value should be above that for which surface coating due to Brownian effects might occur 
(refs. 9 and 10) and below that for which deposition due to overloading and inertial effects 
might result (refs. 11 and 12). 
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